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ABSTRACT 

A Monte Carlo analysis was conducted to determine a theoretical 
Apollo Saturn V launch release design bending moment at a critical 
location on the spacecraft structure. Inputs to the analysis were de- 
fined by their distribution functions and selected by a random number 
generator. Three-dimensional structural dynamic characteristics of the 
Apollo Saturn V vehicle were used' in a forced response program to deter- 
mine the spacecraft loads for 200 cases. The results indicate that the 
maximum bending moment is lognormally distributed. In addition, the 
design bending moment was compared to one obtained using the current 
design method of root— sum-squaring the results of several separate cases . 

' * v I • I , 

This investigation shows that a Monte Carlo approach to determine 
a lift-off design load is feasible, and that the result obtained through 
the current design method is conservative. 
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CHAPTER I 


INTRODUCTION 

Although, in appearance the lift-off of a giant space vehicle is 
slow and majestic, structurally it can he one of the most violent and 
critical phases of the vehicle flight. Prior to lift-off, steady winds 
and gusts create thousands of pounds of lateral drag. Von Karman vortices 
(vortex shedding) occur at critical windspeeds and cause the vehicle to 
sway in a direction normal to the wind. The unsymmetric buildup of engine 
thrusts causes large lateral forces and torques, as well as longitudinal 
forces, to he applied. The application of these wind and thrust forces 

I 

results in static and dynamic deflections of the vehicle and large loads 
in the hold-down structure. The quick, release of these constraining 
loads induces iarge lateral and longitudinal structural oscillations in 
the structure. The design loads on a spacecraft due to these oscillations 
are an important design consideration and are difficult to establish 
because of the many variations in sequence times, thrust buildup rates, 
and maximum thrust. Although considerable statistical information about 
these parameters exists, a meaningful application to determine design 
loads has not yet been established. For example, the design lateral loads 
for the spacecraft are now determined in the following manner (Refer- 
ence l). First, a design windspeed versus altitude profile and a corre- 
sponding vortex shedding load are selected, and the bending moment 
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resulting from a lift-off in that •wind is calculated. Next, a combina- 
tion of thrust buildup forces is selected which produces maxi mu m exci- 
tation of the cantilevered vehicle. The bending moment resulting from 
lift-off during this excitation is calculated. Finally, the results of 
the two cases are root-sum-squared to determine a design bending moment. 

This design method has several drawbacks. First, design loads 
evolve from a selection of worst— worst inputs . Because of the complex 
nature of the problem it is not possible to determine whether or not this 
is a conservative approach. Second, the design case is hypothetical 
since it is not based on a realistic combination of inputs. Therefore, 
the probability of its occurrence is unknown. These drawbacks plus the 
availability of statistical input data, have led to the evolution of the 
Monte Carlo approach of calculating lift-off loads . 

The objective of this investigation is to develop a means of 
statistically analyzing spacecraft loads at lift-off through the use of 
the Monte Carlo method. This method is a solution of a problem in which 
the input data for many cases are selected at random through the use of 
a table of random numbers and probability distributions of the inputs. 
Selected outputs of the problem are then analyzed statistically by con- 
structing their probability distributions from the results of many cases. 
This investigation will apply the technique to determine the Apollo 
Saturn V vehicle spacecraft loads at thrust buildup and launch release. 
Probability distributions of windspeed and direction , engine start times , 
buildup' rates , maximum thrust levels and engine alignments, and launch 
release time will first he established. Using a table of random numbers, 
data for 200 cases will be selected from the probability distributions. 
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The data for each case will he used to set up force time histories and 
sequences of events. Using a digital computer solution, these forces 
will then he applied to a structural mathematical model of the Apollo 
Saturn V and maximum spacecraft structural loads calculated for each 
case-. The peak loads for the 200 cases will then he ordered from low- 
est- to highest, a mean and standard deviation established, and a prob- 
ability distribution plotted. • This distribution will -then be used to 
establish a theoretical design value which will he compared both in 
magnitude and meaning to one obtained by the current design method of 
solution. 



CHAPTER II 


EQUATIONS OF MOTION 

The method used to solve for the dynamic response of the space 
vehicle utilizes the standard normal mode equations with viscous, damping. 
A basic description of the method follows, and a more rigorous derivation 
may he found in Reference 2. 

Consider the Apollo Saturn V space vehicle to he an n degree-of— 
freedom lumped parameter system with mass matrix [M] , stiffness ma- 
trix [k] , damping matrix [C], and the column matrix of external 
forces {F(t)} all expressed in the w coordinate system. The differ- 
ential equations of motion in the w coordinate system take the form 

+ [C]{w} + [K]{w> = {F} (l) 

In the normal mode method, mode shapes and frequencies for a system 
defined hy Equation (l) are found from the matrix equation 

{ [K] - [AJ }{w> = {0} (2) 

The solution to Equation (2) in terms of mode shapes and frequencies is 
found by making the matrix transformation 

{w} = [<£]{£) (3) 

in which each column of [$] is a modal column of the system and {?} 



5 


are normal coordinates. Differentiating Equation C3) with, respect to 
time leads to 

{w} = UHf} ' Ik) 


and 

’ {w> = Mm (5) 

Substituting Equations (3), (4) , and (5) into Equation (l) and premulti- 
plying by the transpose of Equation Cl) becomes • * 

ufMUHu + [4»] T icir<i>]{|} + E^uawm = c<i>] T {F> to 


But the orthogonality relationship among normal modes is expressed by 


ihr cM M*ji = ° 

i # 0 

(T) 

Itf »{+;, > - 0 

i T 5 <5 

(8) 


Therefore 

[<t>J T [M]{(j.] ={M g ] (9) 


where 


[M g ] 


is the matrix of generalized masses, and 


* [“r 2 ][“g] = [ K g] 


( 10 ) 
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Comparison, of the triple matrix product [<j>] [C][<j>] with Equations (9) 
and (10) shows that the product results in a diagonal matrix only when 
[ C ] is proportional to either. [M ] or [K] . Making the assumption 
[C] = then 

E^flCJW * 2?[M g ] (11) 

Expressing [<j>] {P} by |F g j., the generalized forces, and substituting 
Equations (9), (lO), and (ll). Equation (6) becomes 

[ M g ]{«r} + Z i M g ]!«ri + [“r 2 K]) 5 rl * ) F g} Cl2) 

Equation (12) then is the standard form of the normal mode equa- 
tions -which can be used to calculate the dynamic response of the space 
vehicle. A description of the structural model is included in Chap- 


ter III. 



CHAPTER III 


STRUCTURAL MODEL 

The structural model used in the analysis is 'a l80-degree-of- 
freedom lumped mass model of the Apollo Saturn V vehicle (Figure l). 

Fifty lumped masses 'were configured to match the mass properties of the 
fourth Apollo Saturn V vehicle (Apollo 9 ) • A composite total vehicle 
lumped mass model was constructed by North American Rockwell' Corporation 
from mass and stiffness matrices of the spacecraft provided by NASA 
Manned Spacecraft Center and of the launch vehicle provided by NASA 
Marshall Space Flight Center* 

The lumped mass models of each of the components making up the 
entire vehicle are shown in Figures 2 through 6. For reference purposes, 
the 50 lumped masses (hereafter called nodes ) are numbered consecutively 
from the top of the spacecraft to the base of the launch vehicle. 

Figure 2 shows the launch escape system (LES) and command module (CM) 
portions of the spacecraft. The service module (SM) , housing spacecraft 
life support systems,- spacecraft propulsion system, and fuel and oxidizer 
tanks, is illustrated in Figure 3. Figures U and 5 are the lumped mass 
descriptions of the lunar module (LM) and its housing and the spacecraft 
IM adapter (SLA), respectively. The nodal description of the three-stage 
launch vehicle is illustrated in Figure 6. The coordinates' of the entire 
lumped mass system are listed in Table I and are consistent with the 
coordinate system of Figure 1. Also shown in the table are the degrees 
of freedom retained for each mass point. 
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Orthogonal modes and frequencies for the launch vehicle were cal- 
culated "by Worth American Rockwell Corporation using the Givens Method 
(Reference 3 ). Two sets of modes were calculated and subsequently used 
in the analysis. The first set is for the ignition phase, wherein the 
vehicle is still held down on the launch pad (cantilevered). A total of 
91 cantilever modes were calculated, 30 of which are used in this study. 
The second set of modes calculated were for the lift-off or free— free 
configuration. A total of 84 free-free modes were calculated, 30 of 
which are used in the study. Modal damping equal to 1 percent of criti- 
cal damping is used in the analysis for all cantilever and free-free 
modes . 



CHAPTER IV 


LOADS EQUATIONS 

The Apollo spacecraft is a very complex structure requiring 
complicated stress analyses to trace load paths for a given set of body 
loads. As indicated in the introduction, however, the objective 1 >of this 
investigation is to derive a method for determining a realistic set of 
basic body loads for the design condition of lift-off. A description of 
the stress analysis stemming from the calculation of body loads is there- 
fore not within the scope of this work. 

Hereafter, the term "body loads" will refer to the combination of 
shears , axial forces , and bending moments existing at a particular loca- 
tion on the spacecraft and resulting from externally applied forces and 
inertia forces . Although a running distribution of body loads is neces- 
sary to structurally design the spacecraft , such a distribution usually 
evolves from the calculation of body loads . The critical locations are 
generally the interfaces where major portions of the spacecraft are mated 
Examples on the Apollo spacecraft are the LES/CM, CM/SM, SM/SLA, and the 
IM/SLA interfaces. Although it is necessary to establish body loads at 
all of these interfaces , this investigation will consider only the CM/SM 
interface as a typical and critical example. Figure 7 illustrates the 
math model lumped mass representation of the spacecraft above the CM/SM 
interface. The body loads to be calculated are those acting upon the SM 
side of the interface and their signs are considered positive in the 


direction shown. The loads are calculated from a summation of externally 
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applied forces and inertia forces at>ove the interface. Externally ap- 
plied forces above this interface consist entirely of wind-induced lateral 
forces. The inertia forces result from rigid and elastic body accelera- 
tions discussed in the previous chapter. Note that rotational degrees 
of freedom were not included in the modeling of the LES and are therefore 
not included in inertia terms of the load equations. The body loads at 
CM/SM interface are claculated from the following equations. 

I. AXIAL FORCE 


S 


x 



( 13 )' 


i 

where the m. f s are the lumped masses and the x. 's are their total 
1 i 

rigid and elastic body acceleration components in the x direction. 

A positive force indicates a tension axial -force at the interface. 


II. SHEAR FORCES 


S 

y 




m i y i 


i=l 


S 

z 


-2 



i-A 


i=l 


( 1 *) 


( 15 ) 
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■(16) 


where F. and F: are the y and z components, respectively, of 

V 1 z 

externally applied forces and S„ is the total shear at the CM/SM inter- 

X\ 

face. 


III. MOMENTS 


. .. n*i 


M 


x 


- Z 3 S + B 3 Z #3 - V3 2 3 

X X 


(IT) 


where M is the torsional load about the x axis at the CM/SM inter- 
x 

face, I_ is the mass moment of inertia of the command module, and 
x 

0 is its rotational acceleration about the x axis . The second and 
J x 

third terms in the equation result from the slight offset of the CM cen- 
ter of gravity from the x axis. There are no externally applied tor- 
sional forces at lift-off. 


M 

y 


T>i ( x i - X CM/SM) - J 3 S 3 + 2'Vi( x i " X CM/S M ) - "3^3 
i?l Z y 7 i=l 


( 18 ) 
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x CM/Sm) + m 3 y 3 X 3 



moments, respectively, at the CM/SM interface. The term (x^ - 
represents the x-axis distance from each lumped mass to the interface. 
The final terms in Equations (l8) and ( 19) again result from the offset 
CM mass. 



CHAPTER V 


INPUT DATA 

As noted in the introduction, the purpose of this thesis is to de- 
termine by means of a Monte Carlo analysis the spacecraft structural 
loads resulting from 200 separate lift-off cases. Each case is composed 
of a set of discrete force time histories applied to the structural model 
and the model response to the forces. This section ■will describe how 
discrete forces for each case are developed from both fixed data and 
those data described by probability distributions. The data may gener- 
ally be classified as sequencing, environmental, or propulsion. 

I. SEQUENCING DATA 

An important aspect of the launch release problem is the sequence 
of events. Variations in ignition or launch release times can have con- 
siderable impact on the ensuing structural loads. 

The nominal lift-off sequence for the Apollo Saturn V vehicle is 
shown in Figure 8. The sequence begins at -5.U seconds by an accurate 
timing device, although nothing significant to the structural simulation 
occurs until 0.0 second. Because the Saturn V first stage has five en- 
gines, the ignition signals are planned nominally to have as little dy- 
namic effect as possible on the launch release. This is accomplished by 
staggering the start times to minimize longitudinal excitation and 
torques or side forces applied by the thrust buildup. Figures 8 and 9 
illustrate how this is accomplished. Engine 5, the center engine. 



receives the ignition signal at O.k'J second. The remaining engines are 
nominally ignited in diametrically- opposing pairs at 0.3-second intervals 
thereafter. At 3.5 seconds the chamber pressures of all five engines are 
sensed, and if all have attained 90 percent of nominal thrust, a launch 
commit signal is given to retract the hold-down arms constraining the 
vehicle to the pad. These arms nominally release 0.22 second after the 
commit signal, or at 3.72 seconds. 

Although the timer is exceptionally accurate , the sequence is 
never exactly nominal due to time variances in actual engine ignitions, 
and hold-down arm retraction after the signal has been transmitted'. 

Enough data, however, have been attained through tests to describe the 
probability distributions of these time variations. Variation in igni- 
tion time of each engine is described by its probability distribution in 
Figure 10. The distribution is normal with a mean of 0.0-second time 
variation from the nominal sequence time and a standard deviation of 
0.055 second from the mean. The probability distribution of launch 
release delay (time from launch commit to hold-down arm retraction) is 
shown in Figure 11. As indicated, the nominal launch release time is 
3.72 seconds, which represents a mean of 0.22 second in launch release 
delay. 

XI. PROPULSION DATA 

A means of statistically describing the thrust buildup of the 
Saturn V engines at first appears remote because of variation in igni- 
tion time, buildup slopes, and maximum thrust. However, a close 
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investigation of the buildup characteristics reveals a characteristic 
shape to the curve. It is this characteristic or nominal shape 
(Figure 12) and statistically described variations from it that permit 
thrust buildups to be constructed for the Monte Carlo study, 

The buildup curve 1 of Figure 12 is characterized in the following 
manner. The time t^ is determined from the sequence statistics al- 
ready described. The first portion of the buildup (from A to B) varies 
little from engine to engine. The buildup from B to C is marked by a 
rapid buildup rate which varies from engine to engine. Once again 

the buildup is marked by a characteristic plateau (C-D) before entering 
another rapid buildup portion (D-E) with rate • The last portion of 
the thrust buildup (E-F) is a leveling off to the maximum thrust. The 
probability distributions used with ignition sequence distributions to 
construct the thrust buildup curves are' shown in Figures 13, l4, and 15 . 
Figures 13 and l4 show the distributions of the rates R^ and R^ » re- 
spectively, Figure 15 gives the distribution of maximum thrust. 

Another important factor in the construction of thrust magnitudes 
is thrust direction. All five engines are nominally directed so that 
their thrust at ignition will act parallel to the longitudinal axis of 
the vehicle. Through an accurate measurement of many engine alignments 
after vehicles have been stacked, enough data have been gathered to con- 
struct a statistical distribution of engine alignments. These align- 
ments are expressed in terms of two actuator positions for each engine 
(Figure 9), The vectors indicate the direction of a positive engine 
alignment for each actuator position. Lateral components of thrust 



1 6 


vectors would then "be directed in the opposite direction. The available 
data on alignments indicate that they are similar for both actuator - 
positions. The data are represented by its probability distribution in 
Figure 1 6 which is used for both actuator positions and for each of the 
five engines. 


III. WIND DATA 

Probably the most important influence on the launch release lat- 
eral loads is that of surface winds. The lateral drag and wind-induced 
vortex shedding are largely responsible for the lateral loads at the 
base of the vehicle. It is the quick release of these constraining 
forces that causes the dynamic response of the vehicle at lift-off. 

The wind data are based on ill years of anemometer data collected 
in the Cape Kennedy area. In order to establish the wind forces on the 
vehicle, it is first necessary to establish wind direction, gust char- 
acteristics, and a windspeed-altitude profile. The windspeed-altitude 
profile is established by an empirical formulation. Because the earth 
surface exerts a frictional force on the lower layers of the atmosphere, 
the speed-altitude profile is usually determined from the "power law" 
equation 



where Vy is the vindspeed at height h and Vy is the windspeed at 

R 

The exponent P is a function of windspeed and 


a reference height h^. 
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ground roughness. For design purposes and for this thesis a value of 
0.2 is used. An example of this windspeed-altitude profile is shown in 
Figure 17(a). 

In, addition to the steady-state wind profile, a time history gust 
shape must he assumed. This gust shape (Figure 17(h)) is represented by 
a wedge with linear increase to the peak wind in 2 seconds and a linear 
decay to steady-state wind in 2 seconds . Substantial data indicate that 
a design gust factor of 1.4 for all altitudes is advisable and is used 
in this study. Actually, since it is difficult to determine steady winds 
from anemometer data, there are more statistical data for peak winds. 
Therefore, a peak wind at a reference height is first selected from its 
probability distribution. This wind is then divided by 1.4 to determine 
a steady-state wind at the reference height. Using this information with 
Figure 17(a) results in a windspeed-altitude-time profile. The prob- 
ability distribution for peak winds at an altitude of 60 feet occurring 
during a 1-hour exposure period for an annual reference period is shown 
in Figure 18. Although these data provide enough information to es- 
tablish a wind profile, an important effect on the loads is the actual 
time at which the peak wind occurs. After a discussion of the problem 
with the wind analysis group at Marshall Space Flight Center , a decision 
was made to use the nominal lift-off time as the mean peak wind time 
and assume the time to be distributed normally with a standard deviation 
of 0.5 second from the mean. The probability distribution for the time 
of wind peak is shown in Figure 19. The remaining characteristic needed 
to define the profile (wind direction) is described in terms of its 
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probability distribution in Figure 20. Another important effect due 
to the wind is that of Von Karman vortex shedding. This effect results 
from "shedding" of vortices at critical windspeeds . When the frequency 
with which the vortices are shed is close to a natural frequency of the 
cantilevered vehicle, dynamic motion is excited. The extent of vortex 
shedding on the Saturn V has been estimated from wind tunnel tests and 
the dynamic characteristics of the vehicle. The magnitude has been 
expressed in terms of bending moments at the base of the vehicle and is 
a function of both steady-state wind and direction. Because of the 
difficulty in expressing vortex shedding forces, an approximation must 
often be assumed. For this study, vortex shedding is represented by a 
force acting laterally on the vehicle at its center of pressure. Its 
magnitude is that required to produce the test-determined base bending 
moment for a given steady-state wind and direction. Although the re- 
sultant vortex shedding loads is a low-frequency sinusoid (f = 0.3 Hz), 
it is expressed conservatively as a constant force at the center of 
pressure normal to the wind direction at the time of lift-off. The plus 
or minus direction of this normal force is determined randomly by the 
flip of a coin. The vortex shedding force versus steady-state windspeed 
and direction curves are shown in Figure ‘21. 



CHAPTER VI 


METHOD OF SOLUTION 

The determination of launch release structural loads is compli- 
cated not only by complex input data, hut by the dynamic transition of 
the structure from a cantilevered to a free- free state. Although several 
approaches are possible for the solution to such a problem, the following 
was chosen. As described in Chapter III, two sets of orthogonal modes 
and frequencies were first calculated — one set for the vehicle canti- 
levered from the launch pad and the second for the same vehicle free- 
free. As illustrated in Figure 22, the external force time histories are 
first applied to the cantilevered vehicle and the dynamic response and 
loads are recorded until lift-off. The time histories of the constrain- 
ing shear, axial force, and bending moment at the base of 'the launch ve- 
hicle are also calculated during this period. 

The next step (Figure 22) is to apply to the free-free structure 
the same external forces plus the calculated base constraining forces 
(equal in magnitude but opposite in sign). In this way, the cantilever 
response is simulated by the free-free structure and constraining forces. 
The free-free simulation is allowed to continue until after lift-off by 
removing the constraining forces at the time of launch release. The con- 
straints are removed by a linear decay to zero in 0.02 second. This ap- 
proximates the actual time necessary to remove all constraining forces 
from the hold-down arms. The CM/SM interface loads are then calculated 
from the free-free response. 



I. COMPUTER SOLUTION 


Until now the discussion has teen entirely theoretical without re- 
gard to the most practical aspect of the problem, the mechanization re- 
quired to apply the theory and to obtain the desired results. This is 
achieved through the use of three electronic computer programs pro- 
grammed for the Univac 1108 digital computer. Two of the programs were 
developed for general purpose studies by NASA contractors and are used 
without modification in this analysis. The third program was developed 
specifically for this study. 

II. RANDOM INPUT GENERATION 

The random input generator program (Reference 4) accepts as inputs 
the probability distributions described in Chapter V and, by use of a 
random number generator, produces the desired input for each case. The 
program is an operational general purpose program. As such, its oper- 
ation and theory in selecting truly ’’random" numbers are accepted with- 
out detailed discussion. The program is denoted by "RIG." 

III. FORCING FUNCTION GENERATOR 

Developed especially for this study is the "SETUP" program. The 
program accepts as input the data randomly generated by "RIG." It applies 
the theory described in Chapter V to this data and outputs, on cards, 1 
of the force time histories necessary to complete a forced response case. 
The program was checked with desk calculator computation. Detailed de- 
scription of ’’SETUP" is given in Appendix B. 



XV. FORCED RESPONSE PROGRAM 


The third program employed in the analysis is that used to solve 
the equations of motion (Chapter II) for dynamic responses and the loads 
equations (Chapter IV ) for CM/SM interface tending moment « "'The* program 
is called "FLAP" (Flight Loads Analysis Program). This program was de- 

e 

veloped hy North American Rockwell Corporation Ly modifying the program 
"DEMR" (Reference 4) to include the loads equations of ChaptermlVn- The 
program accepts as inputs either the cantilevered or free-free orthogonal 
modes discussed in Chapter III. The modal data input on magnetic tape 
was provided by North American Rockwell. The program also accepts the 
force time histories from "SETUP" on cards. Using these inputs", the 
program calculates rigid and elastic motions and accompanying time his- 
tories of "body loads. 

The flow chart showing the use of the aforementioned computer pro- 
grams is shown in Figure 23. 



CHAPTER VII 


NUMERICAL PROBLEM 

To illustrate best how the force time histories are set up from 
the input data and how the launch vehicle responds to these forces, a nu- 
merical problem has been selected and the details of its solution are 
presented in this chapter. 

The case selected as an example results in a maximum CM/SM bending 
moment of 1.2 x 10^ in-lb,- a load which was exceeded no more than 5 per- 
cent of the time in 200 cases. Table II illustrates the random input 
data selected for this case by "RIG" from the probability distributions. 
The first five numbers are the variations from the nominal ighition times 
selected from Figure 10 for each of the five engines. Engines 2, 3, and 
5 ignited early while engines 1 and 4 ignited late. The next five data 
points selected from Figure 13 are the buildup rates R^ for engines 
1 through 5 , respectively (Figure 12). Note the rapid buildup rate of 
engine number 2. The next five numbers in the table are the buildup 
rates (Figure 12) for each of the five engines which were selected 

from Figure l4. The maximum thrust values for each of the engines con- 
stitute the next five data points in Table II. The distribution from 
which they were selected is shown in Figure .15 . These randomly selected 
ignition times , slopes , and maximum thrust , together with the fixed 
thrust characteristics of Figure 12, result in the thrust buildup curves 
shown in Figure 24. The randomly selected alignment angles of each of 
the engines actuator positions were selected from Figure 1 6 and are also 


shewn in Table II. 
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They are used with the thrust buildup curves and the "SETUP" pro- 
gram to establish propulsion force time histories to be applied to the 
structural model. These forces and moments are applied in "FLAP" to the 
translational and rotational degrees of freedom of the launch vehicle 

i t 4 * 

thrust structure. 

The wind characteristics, the next three data points in Table IX, 
were selected randomly from their probability distributions (Fig- 
ures 18, 19, and 20). The wind at the 60-foot reference height peaks 
at 28. 4 knots from a steady-state condition of 20.2 knots (Chapter V). 
The wind peaks 0.24 second after the nominal launch release time and is 
blowing from an easterly direction (6 = 85.5° from the north). Although 
the wind profile is now defined in a form shown in Figure 17 , the 
lateral drag forces must be distributed along the structural model of 
the vehicle in lumped fashion. Drag loads are applied to l6 locations 
along the vehicle where y and z lateral degrees of freedom have 
been retained in the structural model. "SETUP" uses the wind direction 
speed and time of peak to break the drag into 32 force time histories, 
l6 each in the y and z directions at these nodes. A vortex shed- 
ding lateral force of l4,000 pounds was selected from Figure 21 using 
the peak wind and wind direction data -already selected. As indicated 
in Chapter V, this force is normal to .the drag force and is applied for 
convenience to the structural node closest to the vehicle center of 
pressure. One of two possible normal directions was selected randomly, 
resulting in the vortex shedding components. The only remaining data 
point to be randomly selected is the launch release delay. The delay 
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selected £rom. the distribution (Figure ll) was 0.22 second, resulting 
in launch release at 3. 72 seconds. Together with the modal data de- 
scribed in Chapter III, this completes the required inputs to the "FLAP" 
program. Because forces such as wind drag are not zero at the initi- 
ation of the case, an initial condition subroutine of "FLAP" is employed. 
This subroutine calculates the modal deflection of each mode required 
to balance initial forces such that there are no initial accelerations. 
The vehicle response at lift-off to the force time histories and release 
for this case is shown in Figures 25 and 2 6 . Figure 25 illustrates the 
axial (x) and lateral (y and z) components of acceleration of node 3, 
which is that of the structure of the command module (Figure 2) ; 

Figure 26 shows the time histories of the CM/SM interface axial force and 
pitch and yaw components of bending moment. The response before lift-off 
is small and is not shown in the figures. 



CHAPTER VIII 


RESULTS 


I. MONTE CARLO RESULTS 


The peak resultant "bending moments ‘for 200 cases were first tabu- 
lated and then ordered from lowest to highest as indicated in Table III. 
Also shown in the table is the percentage of cases that were less than 
or equal to each case. Thus, for example, 75 percent of the peak bending 

g 

moments were less than or equal to O.656 x 10 in-lb. From the table a 

■rli 

cumulative distribution function of Mg (maximum bending moment x 10 ) 

max 

was constructed as described in Appendix A. This function is shown in 

Figure 27- The distribution function is plotted on normal-probability 

paper which displays data with a normal distribution as a straight line. 

The closeness of a plot to a straight line on normal-probability paper is 

a measure of the closeness of the distribution of the data to a normal 

distribution. Figure 27 is indicative that the- distribution of Mg 

' g max 

does not closely fit a normal distribution. The x test for good- 
ness of fit (Appendix A) was therefore not attempted with distribution 

of Mg . Instead, the distribution function for log (Mg ) was cal- 
max max 

culated and is shown in Figure 28. Note that log (Mg ) straightened 

max 

the tail of the cumulative distribution function of Mg . The sample 

max 

mean and standard deviation for log (Mg ) were calculated and the 

max 2 

assumed normal distribution function plotted in Figure 28. The x 

test described in Appendix A was applied to investigate the hypothesis 

that log (Mg ) is normally distributed. The results of the test 
max 
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(Appendix C) indicate that the hypothesis of a normal distribution is 
acceptable, The maximum bending moment at the CM/SM interface therefore 

\ 

follows the normal-probability law after the logarithmic transformation. 

Such functions are said to be lognormally distributed. 

Before determining the design load, an acceptable launch success 

criterion must first be established. Assume that the design load shall 

be that which would not be exceeded 99 times out of 100 launches. The 

99-percent value of log (M^ ) from the normal distribution in Fig- 

max g 

ure 28 is 0.4. This represents a bending moment of 1.5 * 10 in-lb. 

£ 

Thus, 1.5 x 10 in- lb is the point estimate for a load which should not 
be exceeded more than 1 percent of the time. The precision of the point 
estimate may be evaluated by means of a confidence limit. The theory 
required to justify the confidence limits for point estimates may be found 
in advanced statistics tests and is beyond the scope of this thesis. 

It is mentioned only to indicate that in practice , a design load should 
be determined, from both the -load distribution function and such a confi- 
dence limit. Since what is sought is a design load, a one-sided confi- 
dence limit would be required. The one-sided confidence limit for a given 
load level determines a load quantity which is almost always greater than 
the given load level. For example, consider the probability statement 

Pr pi 99 £Mj = 0.95 

This one-sided confidence statement means that the true, unknown 
99-percent load is less than or equal to M with 95-percent confidence. 
The procedure for' calculating M with an assumed knowledge of the 
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underlying probability distribution (in this case normal) is called 

parametric. The parametric 95-percent confidence limit for the 99- 

percent point estimate of log ) is 0.516. Taking the antilog, 

max 

M = 1.68. 

^max 

Thus, if this confidence limit were considered, we could say that 

g 

our point estimate of a 99-percent load is 1.5 x 10 in-lb, and in addi- 
tion we have 95-percent confidence that the real 99-percent load is less 
6 

than 1.68 x 10 in- lb. The consideration of confidence limits is common 
in engineering design and therefore a likely choice for the design bend- 
ing moment would be 1.68 x 10^ in- lb. Before comparing this to a design 
load obtained by the current method of solution, a discussion of the ac- 
companying CM/SM interface shear and axial loads is necessary, for al- 
though this interface is bending moment critical, it is impossible to 
design a structural interface without a complete description of loads. 

The results of the 200 individual simulations indicate that the range of 
axial loads for all cases were all similar and varied from 22,000 to 
30,000 pounds compression at the time of maximum moment. Therefore, the 
structure should be designed for a bending moment of 1.68 x 10 in- lb and 
this range of axial force. While the interface structure is sensitive to 
ranges in bending moment and axial load, it is less sensitive to varia- 
tions in lateral shear load. Nevertheless , design shear load must ac- 
company the combination of bending moment and axial load. The CM/SM 
shear load ranged from 0 to 10,000 pounds in the 200 cases. Because of 
the large vehicle bending characteristics , bending moments were often 
accompanied by little or no shear. Therefore, the designer must consider 
the entire range of shears together with the range of axial loads -and the 


design bending moment. 
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Although this thesis considers only the CM/SM interface design 
load, a similar procedure would he followed for all other critical inter- 
faces. 


II. CURRENT DESIGN METHOD RESULTS 

In order to compare the Monte Carlo results with those obtained 
by the current method of solution, a "design” case was run using the same, 
structural model and the current method. A brief description of the 
method follows. 

The design bending moment is obtained by root-sum-squaring the time 
histories of loads from three separate cases. The first case calculates 
launch release loads resulting from application of a design wind profile. 
The second ease is the launch release response of the vehicle to a design 
unsymmetrical thrust buildup. The third case is the response to design 
engine alignment forces. The third case is a second order effect compared 
to the first two and was combined with case two for this thesis. 

The current design criteria specify that the vehicle must be de- 
signed to launch in 95-percent winds. These wind loads are for conserv- 
atism to be accompanied by the maximum vortex shedding load. This, 
resulted in a peak wind of 25.2 knots and a vortex shedding load of 
57 5 700 pounds. For case two, the criteria specify the maximum unsymmet- 
ric thrust buildup to be constructed in such a way that a maximum thrust 
induced torque be applied to the base of the vehicle by a combination of 
the ±99. 73-percent engine start times, buildup slopes, and maximum thrust 
values for the five engines. The design alignment is specified by using 
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the 3cf values of the normally distributed thrust alignment. The fol- 

4 » 
lowing equation is used to define the alignment for each of the five 

engines . 


i W\ 


Design engine deflection = — , 0 




3o 


where n indicates the number of engines (5) and 0^ i- s the 3o 
alignment defined in the input data. Both cases were run using the same 
structural model and programs used for the Monte Carlo cases. The maxi- 
mum moment resulting from the root-sum-square of the two cases was 

g 

2.97 x 10 in-lb. The current method specifies the accompanying design 
shear load be a range from 0 to maximum shear resulting from the design 
case, or 0 to 11,000 pounds. The current method also specifies that the 
range of axial load be defined by the peak-to-peak axial forces at the 
time of maximum moment. This resulted in a range of 22,000 to 
30,000 pounds compression. 


III. MONTE CARLO AND CURRENT CRITERIA RESULTS COMPARISON 


Although the current design criteria consider a 95-percent wind- 
speed profile, the inclusion of maximum vortex shedding load regardless 
of windspeed or direction and maximum applied torques due to unsymmetric 
thrust buildup results in a design bending moment with a very low proba- 

g 

bility of exceedence. For example, the value of 2.97 * 10 in-lb obtained 
by the current design method was not once exceeded in the 200 Monte Carlo 
cases. Comparing log (2.97 ) » or 1.09, with the normal distribution .in 
Figure 28 shows, in fact, that we would not expect this load to be 



exceeded at least 99*99 percent of the time. This is not surprising, 
particularly in view of the unlikely (if not impossible) combination of 
wind and vortex shedding load; the latter is considered currently to be 
maximum regardless of windspeed or direction. 

Although the bending moment varied considerably, the ranges of 
axial load and shear to be considered were similar for both the current 

/ 

and Monte Carlo methods. 

IV. CONCLUSIONS AND RECOMMENDATIONS 

The foregoing discussion results in the following conclusions. 

a. A Monte Carlo method to determine a lift-off design bending 
moment at the CM/SM interface is feasible and results in a lognormal dis- 
tribution. 

b. The current design method apparently results in a conservative 
design bending moment. 

c. The Monte Carlo Method can erase some of this conservatism 
from the designs and still produce a realistic design load with confi- 
dence . 

In view of the above conclusions and the fact that this thesis 
constitutes a feasibility study only, the following recommendations are 
made . 

a. The Monte Carlo Method should be pursued further as a logical 
approach to the determination of lift-off design loads . 

b". The theory and application of confidence limits for point 
estimates should be investigated thoroughly in such an analysis. 
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c. Attempts should he made to determine correlation between input 
parameters and the resulting loads. For example, greater launch day 
operational capability could be provided if it were discovered that high 
loads result only from a combination of specific windspeed and direction. 
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TABLE I.- LUMPED MASS MODEL COORDINATES AND DEGREES OF FREEDOM 


Degrees of freedom retained 



1 * 209.5 

1*052.5 

3798.5' 

371*9.5 

3586.8 

3614.0 

3665.0 

3716.8 

3614.0 

3665.0 

3716.8 

3611*. 0 

3665.4 • 

3716.8 

3614.0 

3665.4 

3691.1 

3596.0 

3397.6 

3307.6 

3316.3 

3316.3 

3316.4 

3316.3 

3341.2 

3258.5 

3101.0 

2832.0 


0.0 

.0 

-.56 

.0 

.035 

48.31 
48.. 31 
48.31 

14.83 

14.83 

14.83 
-48.31 
-48.31 
-48.31 
-14.83 
-l4 . 83 
-14.83 
.0 
.0 
.0 
.0 

-54.0 

.0 

54.0 

.0 

.0 

.0 

.0 


0.0 

.0 

6.64 

.0 

.221 

6.58 

6.58 

6.58 

47.75 

47.75 

47.75 

-6.58 

-6.58 

-6.58. 

-47.75 

-47.75 

-47.75 

.0 

.0 

54.0 

-54.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 


14 15 

20 21 

25 26 

27 28 

29 30 

32 33 

34 35 

36 37 

39 4o 
4l 42 
43 44 

46 47 

48 49 

' 50 51 

53 54 

59 60 | 

62 63 

67 68 

72 73 

77 78 

80 81 

85 86 87 

91 92 93 

97 98 99 

101 102 ; 103 




56 

: 57' 

64 

65 

69 

70 

74 

75 

• 
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TABLE I.- LUMPED MASS MODEL COORDINATES AND DEGREES OF FREEDOM - Concluded 


Node 

Coordinates, in 

Degrees of freedom 

L ■ 

retained 

X 

Y 

z 

X 

Y 

I 

z 

'0 

X 

1 

: 0 
T • 

i 6 

z 

29 


0.0 

0.0 

io4 


— 

— 

■ 

— 

30 

2832.0 

.0 

.0 

105 

— 

— 

— ; 

B 

— 

31 

2646.0 

.0 

.0 

106 

107 

108 

— 

B 

— 

32 

2747.0 

.0 

.0 

" 109 

no 

m 

112 

113 

ll4 

33 

2519 . 0 

.0 

.0 

115 

116 

117 

ns 


w 1 m 

34 

2387.0 

.0 

.0 

119 

120 

121 

122 
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2117.0 

.0 

.0 

123 

124 

125 

126 

— 

— 

36 

1848.0 

.0 

.0 

127 

128 

129 

130 

— 

— 

37 

1848.0 

.0 

.0 

131 

— 

— 


— 

— 

38 

1848.0 

.0 

.0 

132 

— 

— 

— 

— 

— 

39 

1760.0 

.0 

.0 

133 

134 

135 

136 

— 


40 

1664 .' 0 

.0 

.0 

137 

138 

139 

l4o 

l4l 

142 

4l 

1564.0 • 

.0 

.0 

143 

144 i 

1 

145 | 

146 

147 

148 

42 

1401.0 

.0 

.0 

j 

149 

150 

151 

152 

— 

— 

43 

1156.0 

.0 

.0 

153 

154 

155 

156 

— 

— 

44 

912.0 

;o 

1 

.0 | 

157 

158 

159 | 

160 

— 

— 

1 45 

912.0 

.0 

.0 

161 

— 

— 

— 

— 

— 

46 

! 

757.0 

.0 

.0 

162 

163 

164 

165 

— 

— 

47 

602.0 

.0 

.0 

166 

177 

178 

179 

— 

— 

1 48 

365'. 0 

.0 

.0 

170 

171 | 

172 | 

173 

— 

— 

49 

365.0 

.0 

- .0 

174 

— 

— 

— 

— 

— 

50 

1 

100.0 

.0 

.0 

175 

176 j 

177 

178 

179 

180 



















TABLE II. - EXAMPLE RANDOM INPUT DATA 


Variation in engine 1 ignition time 

0.0017^ sec 

Variation in engine 2 ignition time 

-.0185 sec 

Variation in engine 3 ignition time 

-.0191 sec 

Variation in engine 4 ignition time 

.040 sec 

Variation in engine 5 ignition time 

-.0476 sec 
£ 

Engine 1 buildup rate , R^ 

4.00 x 10 ° lb/sec 
c 

Engine 2 buildup rate, R^ 

11.76 x 10 Ib/sec 
£ 

Engine 3 buildup rate, R 1 

' 4.99 X 10 ° lb/sec 

Engine 4 buildup rate, R 1 

7.64 x 10 lb/sec 

Engine 5 buildup rate', R^ 

5.58 x 10 ° lb/sec 

Engine 1 buildup rate-, R 2 

1.69 x 10 lb/sec 

Engine 2 buildup rate , 

1.73 x 10 ° lb/sec 

Engine 3 buildup rate , 

1.49 x 10 Ib/sec 

s; 

Engine 4 buildup rate, R 2 

1.81 x io b lb/ sec 

Engine 5 buildup rate , R 2 

1.32 x 10 ° lb/sec 

Engine 1 maximum thrust 

1.50 x io° lb 

sr 

Engine 2 maximum thrust 

1.51 x 10 lb 

Engine 3 maximum thrust 

1.52 x 10 lb 
£ 

Engine 4 maximum thrust 

1.52 x io D lb 

z' 

Engine 5 maximum thrust 

1.50 x 10 lb * 

Engine 1 actuator 1 misalignment 

0.00835 radian 

Engine 1 actuator 2 misalignment 

0.00820 radian 

Engine 2 actuator 1 misalignment 

0.00139 radian 

Engine 2 actuator 2 misalignment 

0.00468 radian 

Engine 3 actuator 1 misalignment 

0.00518 radian 

Engine 3 actuator 2 misalignment 

0.00736 radian 

Engine 4 actuator 1 misalignment 

0.00467 radian 

Engine 4 actuator 2 misalignment 

0.00451 radian 

Engine 5 actuator 1 misalignment 

O.OO 685 radian 

Engine 5 actuator 2 misalignment 

0.0029 radian 

Peak vindspeed 

28.39 knots 

Wind direction 

85.5 degrees from N 

Wind peak time 

3.96 sec 




Cumulative 
probability , 
percent 

CM/SM 
interface 
bending moment, 

in-lb x 10 

0.5 

0.181 

1.0 

.189 

1.5 

.196 

2.0 

.209 

2.5 

.214 

3.0 

.215 

3.5 

.231 


Cumulative 

probability, 

percent 


20.0 

‘20.5 

21.0 

21.5 

22.0 


in-lb x 10 


0 . 31*9 

. 31*9 

.352 

.353 

.359 


8.0 

8.5 
9.0 

9.5 

10.0 


25.0 

25.5 

26.0 

26.5 

27.0 


10.5 

11.0 

11.5 

12.0 

12.5 


27.5 

28.0 

28.5 

29.0 

29.5 


13.0 

13.5 

14.0 
ll*.5 

15.0 


30.0 

30.5 

31.0 

31.5 

32.0 


15-5 

16.0 

16.5' 

17.0 


32.5 
33.0 

33.5 
31 *. 0 
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TABLE III.- ORDERED 


Cumulative 
probability , 
percent 

CM/SM • 
interface 
bending moment, 

in-lb x 10 6 

34.5 

0.422 

35-0 

. .428 

35.5 

.432 

36.0 

.435 

36.5 

.436 

37.0 

.436 

37.5 

.437 

38.0 

.437 

38.5 

.441 

39.0 

.446 

39.5 

.447 

4o.o 

.447 

40.5 

.448 

4i.0 

.451 

41.5 

.451 

42.0 

.456 

42.5 

.456 

43.0 

-457 

43.5 

.460 

44.0 

.461 

44.5 

.465 

45.0 

.470 

45.5 

.472 

46.0 

.473 

46.5 

.478 

47.0 

.479 ' 

47.5 

.494 

48.0 

.495 

48.5 

.496 

49.0 

.501 

1 49.5 

.502 

50.0 

.505 

50.5 

.508 

51.0 

.512 

1 51-5 

.514 


LOADS - Continued. 


Cumulative 
probability , 
percent 

CM/SM 
interface 
bending moment 

in-lb x 10" ' 6 

52.0 

0.515 

52.5 

.519 

53.0 

.519 

53.5 

.521 

54.0 

.522 

54.5 

.522 

55.0 

.523 

55.5 

.527 

56.0 

.534 

56.5 

.534 

57.0 

.540 

57.5 

.540 

58.0 

.545 

58.5 

.548 

59.0 

1 

.550 

59.5 , 

.551 

60.0 I 

• 551 

60.5 

.556 

61.0 

.558 

61.5 

.560 t 

i 

62.0 

.562 

62.5 

.563 

63.0 

.564 

63.5 

.565 

64.0 

• 570 

64.5 

• 570 

65.0 

• 571 

65.5 , 

• 572 

66.0 

.583 

66.5 

.585 

67.0 

.587 

67.5 

.587 

68.0 ; 

.588 

68.5 • 

.589 

69.0 

.590 
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TABLE IIL- ORDERED LOADS - Concluded 



CM/SM 


CM/SM 

Cumulative 

interface 

Cumulative 

interface 

probability , 

bending moment. 

probability , 

bending moment. 

percent 

in-lb x io -6 

percent 

in-lb x 10" 6 

69.5 

0.594 

85,0 

0.727 

TO. 0 

.603 

85.5 

.741 

70.5 

.608 

86.0 

.766 

71.0 

.611 

86.5 

.771 

71-5 

.612 

87.0 

.787 

72.0 

.614 

87.5 

.789 

72.5 

.614 

88.0 

.799 

73.0 

.630 

88.5 

.801 

73-5 

.637 

89.0 

.805 

74. 0 

.639 

89.5 

.829 

74.5 

.654 

90.0 

.839 

75-0 

. 656 

90.5 

.848 

75-5 

.664 

91.0 

.858 

76.0 

.671 

91.5 

.891 

76.5 

.671 

92.0 

.894 

77.0 

.682' 

92.5 

.922 

77.5 

.686 

93.0 

1.022 

78.0 

.686 

93.5 

1.023 

78.5 

.686 

94.0 

1.053 i 

79-0 

.687 

94.5 

1.117 

79.5 

.697 

95.0 

1.202 

80.0 

.698 

95.5 

1.213- 

80.5 

.699 

96.0 

1.274 

8l.O 

.700 

96.5 

1.283 

81.5 

,706 

97.0 

1.408 

82.0 

.709 

97.5 

1.535 

82.5 

.714 

98.0 

1.554 

83.0 

.718 

98.5 

2.305 

83.5 

.721 

99.0 

2.455 

84.0 

.726 

99.5 

2.536 

84.5 

.727 

100.0 

2.860 
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Figure 1.- Apollo Saturn V lumped mass model. 




Figure 2.- Launch escape system and 
command module lumped mass model. 




Figure 3.- Service module lumped mass model. 
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Figure 8.- Nominal Saturn V lift-off sequence. 
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Figure 9.- Placement of Saturn V engines. 
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Figure l6.- Cumulative probability distribution 
of actuator 1 and 2 alignments. 
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Figure 22 • - Application of cantilever and free-free models to simulate launch release. 
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Figure 23.- Computer solution flow chart. 
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APPENDIX A 


THE MONTE CARLO METHOD 

A Monte Carlo method is generally described as one in which a 
given problem is simulated by- some -suitable random- process;- More 'spe‘- 
cifically, it may be defined- in application to this thesis' as the sta- 
tistical determination of a- lift-off design load where the input data are 
first selected at- random and' the output data (structural loads) are then 
analyzed statistically. Random input data' are selected and loads cal- 
culated for many cases in order to obtain a sample of ' output- data ade- ’ 
quate to define- its distribution. In order to conduct and understand * f 
the results of- such a study, a few principles of probability are neces-' 
sary and 'are derived -below. 

• Let the one-dimensional set’ A- he such that 

f(x)dx = 1 . f(x) > 0 

for all x in A and f(x) has, at most, a finite number of disconti- 
nuities in every finite interval that is. a subset of A. 



If A -is the sample space of the random variable X, and if the 
probability that X . is in A is defined by f f(x)dx, then X is a 

Ja 



continuous-type random .variable wixh a probability density function 
(p.d.f. ) f(x). The cumulative distribution function F is now defined. 

X 

r x 

F = Pr[X 5 x] = / f(x)dx 

Since f(x) >0, F is a nondecreasing function. 

Suppose now that x is one of the inputs to the problem and is 

described by its distribution function F^. Assume also that to solve 

the problem it is desired to make n simulations requiring n random 

variables X, , . . , , X from the distribution function F . To obtain 
I n x 

a random sample of the variable X, it is sufficient to obtain a random 
sample of a variable X, which is uniformly distributed over the interval 
0 to 1, This follows from the fact that the distribution function F 

x 

of the random variable X is a nondecreasing function ranging from 
0 to 1. Consequently, (y), an inverse function, may be defined for 

values of y between 0 and 1. F x _1 (y) is equal to the smallest value 

of x, satisfying the condition that F > y. 

X 

In terms of the inverse function F -1 (y) to the distribution 

X 

function F of the random variable X, the following theorem may be 
stated: 

Let Y^j Y^, ...» Y^ be independent random variables each uni- 

formly distributed over the interval 0 to 1. The random variables defined 


X 


1 



X, 




are then a random sample of the random variable X (Reference 5). 



Tl 


This theorem permits the selection of n random variables whose - 
distribution function is known by selecting n random numbers uniformly 
distributed between 0 and 1. . - 

The same principles apply also to the analysis of output data. 

The primary statistical quantity of interest in. this Monte Carlo analysis, 
is the probability that a given load will not be exceeded. For example, 
it may be desired that the design load be that which would not be ex- 
ceeded 95 percent of the launches. With a finite sample size it is im- 
possible to obtain the exact 95- percent load, but a point estimate of this 
quantity is possible. This point estimate is obtained in the following 

manner. If the load can be assumed to be normally distributed with mean 

2 

y and variance c then the true 95-pereent load is 


F .95 W + K .95 CT 


The corresponding point estimate of this load is given hy 


F« „ = F + K _.S 


.95 


95 


where 



and 
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— 2 

where F is the sample mean and S is the sample variance. The square 
root of the sample variance S is called the standard deviation and If 
is the number of cases in the sample. The value of the constant K 

• 95 

is 1.645 as obtained from a table of normal percentage points. The above 
are unbiased estimates of the mean and standard deviation of the entire 
population. The sample, variance, however, if calculated from 


s2 = IE( F i- p ) 


•would not be an unbiased estimate of the variance 


of the population. The sum of the squares of deviation of the individual 
sample values from the sample mean results in the minimum sum of squares 
of deviation. The sum of squares of deviation from any value other than 
the mean will he larger than from the mean. Since the sample mean may 
not he identical with the population mean, the sum of squares of deviation 
of the individual sample values from the sample mean would be less than 
the sum of the squares of deviation of the individual sample values from 
the population mean. The standard deviation of the sample, computed from 
the sum of the squares of deviation divided by N would therefore be 
smaller than if the sum of squares had been calculated from the true pop- 
ulation mean. To overcome this bias , the variance estimated from the 
sample is obtained by dividing the sum of squares of deviation by If - 1 
instead of by If (Reference 6). 

The statistical estimates described above are based on the as sump- 

I 

tion that the underlying loads distribution is normal. The validity of 
this assumption can he investigated by means of the x goodness of 
fit test. The If loads should first be ordered from minimum to maximum 
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values and a cumulative distribution function plotted. Then the U 
loads are classified to m mutually exclusive intervals. The. .number of 
observations in each interval is n^, i = 1, 2, . .., m. Then the number 
of expected observations (e^, i = 1, 2, . . . , m) , -which according to the 
theoretical distribution should fall into each interval, are calculated. 
The following statistic is then calculated. 



If the observed and theoretical probability distributions are 
2 

identical, the x test statistic is distributed according to the 
2 

X probability law with m - three degrees' of freedom. The degrees of 

2 

freedom which enter into the evaluation of X are the number of inde- 
pendent observations that are available for its calculation (Reference h). 
Thus m intervals provide m - 1 independent intervals. The expected 
observations depend on the two-point estimates of the mean and variance. 
Thus, the number of independent observations are m-l-2=m-3. 

If the observed and theoretical probability distributions are 

o 

identical, the value of the x test statistic will be small. The 

probability of occurrence of any given value of the test statistic may 

2 

be determined from a table of x percentage points. For example, 
if the calculated value of the test statistic for 12 degrees of freedom 
is 18.5, then the corresponding probability is 10 percent that the vari- 
ation of the observed loads from the theoretical distribution would exceed 



the calculated variation, even if the observed and theoretical distribu- 
tions were identical. For individual data quantities in engineering 
practice, the hypothesis of identical distributions is usually rejected 
at the 5-percent level. 



APPENDIX B 


PROGRAM "SETUP" 

The "SETUP" program is the only computer program developed specif- 
ically for this thesis. Its sole purpose is to accept the randomly 
generated data from "RIG" in card format and setup forcing functions 
acceptable to the "FLAP" program. Symbols used in this program are 
listed at the end of this appendix. 


Forces Required 

"SETUP" is used for this study to setup forcing functions for the 
program "FLAP" for each case. To do this it must accept randomly gener- 
ated data from "RIG" and establish discrete force versus time tables. 
Outputs include F_ , F„, , and F^ which are the x, y» and z, com- 
ponents of thrust acting on the corresponding degrees of freedom of the 

thrust structure. F and F are thrust-induced turning moments ap- 

T li T 5 

plied to the y and z rotational degrees of freedom at the same node. 

F_ and F_ are the y and z components of the vortex shedding force 

T 6 T 7 

applied as a constant to the node closest to the launch vehicle center of 

pressure. F and F^ through F and F are alternately the 

i 8 1 9 X 38 i 39 

y and z components of drag force time histories applied to 1 6 nodes 


along the vehicle. 



76 


Equations 

a. Wind profile, W ± versus 


W= V /l.k 
1 peak 


t ¥ ~ 0*® 
1 


w 2 = «i 


t TT = t . - 2.0 
W c peak 


W = V 
3 peak 


t IT * t . 
W^ peak 


\ = w x 


t tT = t , +2.0 
W, peak 


W 5 = W l 


t TT = 10.0 

W 5 


b. Thrust curves, T. . versus t... 


T. n = F 
ll 1 


t._ = 0.0 
il 


T = F 
i2 2 


t. 0 = t. + At. 

i2 l.i 


T - 
i3 

F 3 

*13 = *12 * 

II 

Eh H 

f 4 

*i4 = *13 * 

T = 

i5 

F 5 

*15 ‘ *14 - + 

T. = 
x 6 

f 6 

t . , r ~ t . _ + 
io i5 

II 

^‘ H 

F 7 

*17 ‘ *16 ♦ 

T i8 ” 

F 8 

*18 = *17 + 

T i9 = 

F 

max. * 
i 

*19 = *12 + 

T = 

ilO 

■ F 
max. 

*110 = 10 -° 


(o.QS x 10 6 /R 2i ^ 
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Because "the "FLAP" program requires identical time tables for all 
forces s the above wind profile and thrust curves are interpellated at 
50 time points. The following force outputs are calculated for each of 
those time points. 

c. Force outputs 

The following equations are used to calculate the thrust forces at 
the jth time point. 

\ = S Tj -J 

F T 2 = (“ cos ^ 5 °' 1 [ T ij(^12 “ 6 ll) ~ T 2j( S 21 + ^22) 

+T 3j( g 31 " S 32) + T kj(hl * %2) + T 5j( 0 52 ” %l)] 

9 

\ - ( - 005 [ T u( 6 ll + e is) * T 2i( P a2 - $ 2l) 

' _T 3j( e 31 + e 32 ) + %( S la ~ hz)* T 5.) ( B 51 + 6 5s)] 

\ - + • T 2J ’ *33) 

- 128 -ff r iJ + T 2 j - T 3J ' %) 

The following equations are used to calculate the constant compo- 
nents of vortex shedding force. 

F Tg “( F vort)( sin 9 w) 

\ = ( F vort)( QOS 6 w) 
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The remaining equations in the program are used to calculate the 

y and z components of aerodynamic drag at each of 1 6 vehicle nodes at 
the jth time point 
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A. 

1 


F J 


max. 

x 


T. 

x 


w. . 
1 J 


R li 


R 2i 


t. 

1 


T. . 


t . . 


peak 


w. 

i 


peak 


"SETUP" Symbols 

coefficient to determine ith drag force from windspeed at refer- - 

2 

ence height, lb/knot 
nominal thrust at jth time point, lh 

maximum thrust of ith engine, selected from Figure 15 

ith force time history applied to a vehicle node, lh 

total drag force at ith node and jth time point, lh 

initial buildup rate of ith engine thrust, selected from 
Figure 13, lh/sec 

second buildup rate of ith engine thrust , selected from 
Figure 14, lh/sec 

nominal ignition time of- ith engine, sec 

thrust of ith engine at jth time point, lh 

adjusted jth time point for ith engine thrust, sec 

time point at -which wind speed will peak, selected from 
Figure 19, sec 

time point associated with ith wind speed, sec 

peak wind speed at the 60-foot reference height selected from 


Figure 18, knots 



8o 


&L wind speed afc the 60-foot reference height at ith time point, 

knots 

B . . misalignment of ith engine in actuator position j , selected from 

Figure 1 6 , radians 

it i variation in ith engine ignition time selected from Figure 10, 
sec 

} direction from which wind is blowing, selected from Figure 20, 

deg 


* 

/ 



82 




Intervals 

Probability 

Observed 

frequency 

✓ *- 

^-Expected 
/ frequency 

: ( n i- e i ) 2 

e. 

7 - 

0.0667 

10 

13.333 

0.833 

8 

.0667 

16 

13.333 • 

.535 • 

9 

.0667 

18 

13.333 

1.634 

10 

.0667 

15 

13.333 

.208 . 

11 

.0667 

12 

13.333 

.133 

12 

.0667 

13 

13.333 

.008 

' 13 

.0667 

10 

13.333 

.833 

14 

.0667 

6 

13.333 

4.033 

15 

.0667 

13 

13.333 

.008 

Total 

1.00 

200 

200.00 

2 1 
X = 17.35- 


2 

The value of the x test statistic is therefore 17.35... Having 

chosen 15 intervals, the validity of the lognormal distribution assump- 

■ >- 

2 w' 

tion must he determined from the x table of • probability levels for 

15 - 3 = 12 degrees of freedom. From the table, the probability of 
2 

finding a x statistic greater than 17.35 for 12 degrees of freedom is 

approximately 16 percent. This means that even if the hypothesis of 

being lognormally distributed were correct , the probability of 
max 2 

finding a x statistic larger than 17*35 is l 6 percent. Since the 
5 -percent level test was easily satisfied, the hypothesis is accepted. 
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APPENDIX C 


THE x 2 TEST OF LOAD DISTRIBUTION 

p 

The x goodness of fit test described in Appendix A was applied 

4 H 4 

to the observed and theoretical distributions of log (M^ ) shown in 

max 

- Figure 28. Fifteen discrete intervals of equal probability were used to 

2 

calculate the following x statistic. 



i=l 


where is the observed frequency (observed number of points to fall 

within the ith interval) and is the expected frequency (expected 

number to fall within the ith interval, i.e. , —r )• 

2 

The following table shows the calculation of the x statistic. 
The observed frequency was taken from the distribution function' of 
Figure 28. 


Intervals 

Probability 

Observed 

frequency 

'Expected 

frequency 

1 

■rag 

1 

0.0667 

11 

13.333 

0 . 4 o 8 

2 

.0667 

9 

13.333 

i.4o8 ' 

3 

.0667 

/ Ik 

' 13.333 

.033 

4 

.0667 

3 13 

13.333 

.008 

5 

.0667 

22 

13.333 

5.634 

6 

.0667 

18 

13.333 ' 

• 

1.634 
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